Multiphoton excitations and inverse population in a system of two flux qubits 
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We study the multiphoton spectroscopy of artificial solid-state four-level quantum system. This system is 
formed by two coupled superconducting flux qubits. When multiple driving frequency of the applied mi- 
crowaves matches the energy difference between any two levels, the transition to the upper level is induced. 
We demonstrate two types of the multi-photon transitions: direct transitions between two levels and ladder-type 
transitions via an intermediate level. Our calculations show, that for the latter transitions, in particular, the in- 
verse population of the excited state with respect to the ground one is realized. These processes can be useful 
for the control of the level population for the multilevel scalable quantum systems. 
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PACS numbers: 85.25.Cp, 85.25.Dq, 84.37.+q, 03.67.Lx 

Manipulation and measurement of the state of scalable mul- 
tilevel quantum systems is one of the key issues for their 
implementation for quantum information processing devices. 
One of the possible realizations of such type of devices is 
based on superconducting Josephson junctions, where the lev- 
els populations can be controlled by the external microwaves. 
When the single or multiple driving frequency matches the 
energy difference between energy levels, the transition to the 
upper level is induced. This provides the tool for manip- 
ulation and characterization of the quantum system - sin- 
gle/multiphoton spectroscopy. 

Recently the spectroscopy of the energy levels in sin- 
gle Josephson-j unction qubits was demonstrated with multi- 
photon excitations by several groups (see for example |H Q, 
La L4, Dll). Moreover, similar multi-photon effects were stud- 
ied in a single superconducting circuit when it behaves as a 
multi-level system JUS El- On the other hand, the multi-level 
quantum system can be realized in coupled qubits, where the 
one-photon driving was studied 011 [IS [II EBB- Such 
two-qubit devices are of practical importance e.g. for building 
the so-called universal set of gates (l5 , Tfjl J3l which is essen- 



tial part for realization of quantum computing 118, JjJ 2(ill . 

In this paper we report the observation of the multi-photon 
resonances in a two-qubit system. The levels population in 
this effectively four-level system can be controlled by exter- 
nal driving field. We demonstrate two types of the resonant 
excitations: direct (from one level to another, when two levels 
are relevant) and ladder-type (via an intermediate level, when 
three levels are relevant). Moreover, our calculations show 
that the inverse population takes place in our driven multi- 
level system. These results are important for new, four-level 
quantum elements which are difficult to realize in quantum 
optics. For example, they can generate extremely large opti- 
cal nonlinearities at microwave frequencies, with no associ- 
ated absorption 12111 . 

The investigated system consists of two coupled supercon- 
ducting flux qubits. A flux qubit is a superconducting ring 



with three Josephson junctions 12211 . For controlling the state 
of the system, the microwave magnetic flux is applied. For 
reading-out the state, the qubits are weakly coupled to the res- 
onant tank circuit. 

The coupled flux qubits are described by the effective 
Hamiltonian in terms of the Pauli matrices: 
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The tunnelling amplitudes Aj and the coupling parameter J 
both are constants in this model. The biases 



i (t) = 2/W$ /W(t) 



(2) 



are controlled by the dimensionless magnetic fluxes / W (t) 



*<(f)/*o 



1/2 through i-th qubit: 

/ (i) (*) = /i + /acSinurf. 



frf 



(3) 



(i) 

Here Jp is the persistent current in the i-th qubit /j = 

$ dc/$o - 1/2, Ac = $ac/$o, and / rf stand for the flux 
introduced by the tank coil. The parameters for the qubit sys- 
tem were obtained from the ground-state measurements [4]: 



A a = 15.8, A b = 3.5, 4 a) $ = 375 ' / p U;$0 = 700 
J = 3.8 in units of h-GHz. 
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ELj/q b ; is the total mag- 
.)/93>j is the symmetrized 



It was shown 123L |24L 12511 . that the response of the mea- 
suring device, the tank circuit, is defined by the magnetic 
susceptibility or the effective inductance of the qubit system. 
Namely, both the phase shift and the voltage amplitude offset 
in the limit of small bias current (/ r f — ► 0) are proportional 
to A = d$ tot /d$ x , where $ tot = 
netic flux and d(...)/d$ x = £<9(. 

(i) 

derivative. Here Li and iv are the geometrical inductance 
and the expectation value of the current in i-th qubit. The 
value Cl 1 = dl$/d<f> x is the inverse effective inductance of 
the i-th qubit. Thus, the tank voltage amplitude is related to 
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When a qubit is in the ground state, its current has a definite 
direction and is defined by the energy derivative with respect 
to dc flux B23I1 . Then the tank voltage amplitude (see Eq. ©) 
as a function of dc fluxes is defined by the second derivative 
of the energy, i.e. by the curvature of the ground-state en- 
ergy. Thus, one expects a dip in the flux dependence of the 
tank res pon se in the vicinity of the energy avoided crossing 
lEM l27ll . In the driven qubits the current has the proba- 
bilistic character. As the consequence, equation includes 



not only the terms with energy derivative, but also the terms 
with the probability derivative J2, H^. When the system is 
driven, its upper level occupation probability is resonantly in- 
creased. The respective derivative is displayed as the alter- 
ation of peak and dip centered around the resonance 1271 12811 . 
Correspondingly, if the tank voltage is plotted as a function 
of two parameters (two partial dc flux biases), the resonances 
appear as the alteration of ridges and troughs. 




FIG. 1: (Color online). Multi-photon resonances in the system of coupled flux qubits. The tank voltage amplitude (a, c, e) and the energy 
contour lines (b, d, f) versus the partial bias fluxes / a and /b for different driving frequencies (from top to bottom): uj/2-k — 17.6, 7.0, 4.1 
GHz. Numbers k—j next to the lines mean that the line relates to the energy difference Ej — E^. The color inset in (f) is plotted for A, Eq. l(4j, 
by solving the Bloch-Redfield equation for the reduced density matrix. 
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In the left panel in Fig. Q] we present the experimental re- 
sults (the voltage amplitude of the tank as a function of qubit 
biases) for the system of two coupled flux qubits. The driving 



frequency lu/2it is 17.6 for (a), 7.0 GHz for (c), and 4.1 for 
(e). The system can be resonantly excited from the level k to 
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the level j when the energy of n photons matches them: 

A^ fc j(/a, /b) = Ej - E k = n ■ hw. 



(5) 



Then along the contour lines defined with this relation the 
resonant structure appears. (Besides, the trough due to the 
ground state curvature is visible at the center, around /b close 
to 0.) The resonances are visualized with the ridge-trough 
line. However, the ridge-trough line is disturbed with increas- 
ing or decreasing the signal; some of these changes are shown 
with white circles. This means changing the effective Joseph- 
son inductance in these points. We argue that this happens 
because of the ladder-type multi-photon excitations to higher 
levels (see also Ref. [29]). 

To understand the experimental results, we plot the energy 
contour lines in Fig.[TJb,d,f) for three frequencies ui, by mak- 
ing use of relation ©. Consider first Fig. Htb)- The black 
and red lines show the positions of the expected resonant ex- 
citations from the ground state to the first and to the second 
excited states respectively. The blue and orange lines are the 
contour lines for the possible excitations from the first and 
from the second excited state to the third excited state. The 
one-photon resonances along the black and red curves are 
clearly visible in Fig.QJa); this was used for the spectroscopi- 
cal study of the system [4]. For better understanding we calcu- 
late the energy levels, by diagonalizing Hamiltonian ([T), and 
plot them at fixed value of the bias flux through qubit a, / a , as 
a function of the bias flux through qubit b, /b, in Fig.|2a); the 
arrows of the length oj/2it = 17.6 GHz and 7.0 GHz (orange) 
are introduced to match the energy levels. The black and red 
arrows in both Fig. |2J a ) an d Fig- Hfb) show the position of 
one-photon transitions to the first and the second excited lev- 
els. The double green and blue arrows in Fig.|2]show the posi- 
tion of the two-photon processes, where the excitation by the 
first photon creates the population of the first and the second 
levels and the second photon excites the system to the upper 
level. We emphasize that these two-photon excitations hap- 
pen via intermediate levels. The position of these expected 
resonances is shown in Fig. 02b) with the blue and green thick 
points. Indeed, there is the change of the signal in Fig. [Tfa) 
in these points. (The two-photon resonance — 2 — 3, shown 
with the blue point, was shown with the numerical simulation 
in Ref. B27IU Moreover, the orange point in Fig. Q2b) stand 

for the ladder-type three-photon excitation, 0—1 3 (with 

one photon to the first excited level and then with two photons 
to the upper level), which is also visible in Fig. [If a). 

In Fig. Qlc) we can see the ridge-trough resonances for the 
driving frequency lu/2it = 7 GHz. Comparing with the con- 
tour lines in Fig. 02d) we easily notice that the resonances 
are one- and two-photon resonant excitations to the first ex- 
cited level. Note that the two-photon resonant excitation is 
direct and happen without any intermediate level in contrast 
to the above described resonances. The higher level excita- 
tions via the first excited state appear due to three- and four- 
photon excitations, as shown with orange (0 1 — 2) and 

pink (0 1 2) dots in Fig.[Tfd). These resonances are 

visible in Fig.QJc). 
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FIG. 2: (Color online). Ladder-type transitions and the inverse popu- 
lation, (a) Energy levels of two coupled qubits at / a = 0.015 versus 
fh (i.e. along the dotted line in Fig.Q](b,d)). (b) Transition matrix 
elements T nm between the eigenstates \E m ) and \E n ). (c) The en- 
ergy levels occupation probabilities Pi under driving with frequency 
u}/2n = 17.6 GHz. 



In Fig. [T]e) the response of the two-qubit system is shown 
for cu/2ir = 4.1 GHz. The lines along the horizontal axis are 
due to direct 1-, 2-, 3-, and 4-photon excitations to the first ex- 
cited state; cf. black lines in Fig.QJf)- Numerous upper level 
excitations via the first excited level appear as the amplifica- 
tion and lowering the signal along these lines. For illustration 
in Fig. Q2 f) beside transitions to the first excited level (black 
line) we plot the red and violet lines which match the ground 
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state and the second and third excited states. 

The transition probability from the state \E m ) into the state 
\E n ) is defined by the absolute value of the matrix element of 
the perturbation. From Eqs. (Q]|3]l these are given by: 



\(E n \v\E m )\ 2 , 
r (b) V P z p z J ' 



(6) 



where we have divided the perturbation by the factor 7p b ^ $o ' 
/ ac . We plot the transition matrix elements in Fig.[2|b) to de- 
scribe the ladder- type excitations shown in Fig.QJb) with the 
green and blue points. The position of the respective reso- 
nances is shown in Fig. |2J a ) an d (b) respectively with green 
and blue arrows and dotted lines. In the left half of Fig. |2|b) 
we plot the transition elements around the transition between 
the three levels Eq, E\, and £3. In the right half we plot the 
transition elements related to three levels Eq, E2, and E3. In 
the latter case the transition element between the higher two 
levels (E2 and E3) is smaller than between the lower two lev- 
els (£0 and E2Y T23 -C 7o2- In contrast, in the former case 
the transition element between the higher two levels (E% and 
E3) is significantly larger than between the lower two levels 
(Eq and E{): T13 3> Tbi- Important to note that in both cases 
the transition element from the ground state to the highest ex- 
cited level, To3> is ver y small. This means that the probability 
of the direct excitation to the highest level is very small - the 
transition is induced exceptionally due to ladder-type mecha- 
nism. 

And finally we calculate the energy level occupation prob- 
abilities under driving. This is done by means of the numeri- 
cal solution of the Bloch-Redfield equation for the qubit sys- 
tem density matrix p l3oTl . The impact of the dissipative en- 
vironment on the qubit system is described by the Redfield 
relaxation tensor. Conveniently the environment is modelled 
as the harmonic oscillator common bath with the ohmic spec- 
tral densities 13011 . Then the coupling of the qubit system to 
the environment is characterized with only one phenomeno- 
logical parameter a, which describes the strength of the dissi- 
pative effects. Important to note that this model accurately de- 
scribes the relaxation between different levels of the systems 
(in contrast to the often used approach with equal relaxation 
rates for all energy levels), since respective relaxation rates 
are essentially different. Then the fitting of the experimental 
graphs is done with Eq. ©, where the expectation value of 
the current in i-th qubit is calculated with the reduced density 



matrix: IK* = — Ip' SpQxr)?'). The result of the calculation 
is presented as the inset in Fig.QJf)- Such fitting gives us the 
parameter for the strength of dissipation a = 0. 1 and the driv- 
ing amplitudes / ac x 10 3 = 4, 6, and 8 for Fig.QJa), (c), and 
(e) respectively. 

The numerically calculated energy level occupation prob- 
abilities are plotted in Fig. |2fc) for oj/2ir = 17.6 GHz, 
/ac = 4 x 10~ 3 , / a = 15 x 10~ 3 versus / b , that is along 
the dotted line in Fig. [TJb). This graph demonstrates two in- 
teresting phenomena, similar to those which exhibit atoms in 



the laser field 13111 . First, the ladder-type resonant excitation 
takes place to the left, where the upper level occupation prob- 
ability P 3 is of the same order as the intermediate level oc- 
cupation probability Pi . This corresponds to the green point 
in Fig. [TJb)- Second, the inverse population appears to the 
right (corresponds to the blue point in Fig. [TJb)). This means 
that the upper level occupation probability Pi is larger than 
the ground state one Po J3, IE l33ll . In our case the inverse 
population at the first excited level is accumulated after relax- 
ation from the third excited level, since relaxation from the 
third level to the first one (shown with solid curved arrow in 
Fig.[2ja)) is larger than both relaxation from the third to sec- 
ond and from the first to the ground state (shown with dashed 
curved arrow in Fig.[2la)). 

In conclusion, the multi-photon resonances in the four-level 
(two flux qubit) system were observed. The multi-photon res- 
onances were of two kinds: direct and via intermediate level 
(when three levels are relevant for the process). Our calcu- 
lations show that in the latter case the stationary state of the 
system with different relaxation rates exhibits the inverse pop- 
ulation. Such effects are relevant for both multi-photon spec- 
troscopy of the system and for the study of new effects in ar- 
tificial multi-level structures. 
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